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Abstract. An exact nonlinear scaling transformation is presented for the local three-dimensional dy- 
namical equations of motion for differentially rotating disks. The result is relevant to arguments that 
have been put forth claiming that numerical simulations lack the necessary numerical resolution to re- 
solve nonlinear instabilities that are supposedly present. We show here that any time dependent velocity 
field satisfying the local equations of motion and existing on small length scales, has an exact rescaled 
counterpart that exists on arbitrarily larger scales as well. Large scale flows serve as a microscope to 
view small scale behavior. The absence of any large scale instabilities in local numerical simulations of 
Keplerian disks suggests that the equations in this form have no instabilities at any scale, and that finite 
Reynolds number suppression is not the reason for the exhibited stable behavior. While this argument 
does not rule out the possibility of global hydrodynamical instability, it does imply that differential 
rotation per se is not unstable in a manner analogous to shear layers or high Reynolds number Poiseuille 
flow. Analogies between the stability behavior of accretion disks and these flows are specious. 
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^ 1. Introduction 



Disks with Keplerian rotation profiles are linearly sta- 
ble by the Rayleigh criterion of outwardly increasing 
specific angular momentum, but are extremely sensi- 
tive to the presence of magnetic fields. A weakly mag- 
netized disk is linearly unstable if its angular veloc- 
ity decreases outward, a condition met by Keplerian 
and almost all other astrophysical rotation profiles 
(Balbus & Hawley 1991). The underlying physics be- 
hind this magnetorotational instability (MRI) is well- 
understood, and the breakdown of the fiow into fully 
developed turbulence has been convincingly demon- 
strated in a large series of numerical simulations 
(Balbus 2003 for a review) . 

Not all astrophysical disks need have the requisite 
minimum ionization level to sustain magnetic cou- 
pling, however. Protostellar disks, for example, may 
have an extended "dead zone" near the midplane on 



Send offprint requests to: 
sbOvirginia . edu 



Steven A. Balbus, 



radial scales from ~ 0.1 to ~ 10 AU (Gammie 1996; 
Fromang, Terquem, & Balbus 2001). This, along with 
other similar cases (e.g. CV disks, cf. Gammie & 
Menou 1998), has led to speculation that there are 
also hydrodynamical mechanisms by which Keplerian 
flow is destabilized (Gammie 1996). 

Before the advent of the MRI, such reasoning 
was orthodox. The pioneering work of Shakura & 
Sunyaev (1973), for example, invoked nonlinear, high 
Reynolds number shear instabilities as a likely desta- 
bilizing mechanism that would lead to turbulence (see 
also Crawford & Kraft 1956). Since, for nonaxisym- 
metric disturbances, there is still no proof either of 
linear or nonlinear stability, this mechanism contin- 
ues to attract adherents (DubruUe 1993, Richard & 
Zahn 1999, Richard 2003). 

Theoretical analysis may have hit an impasse, but 
the intervening years have in fact seen a stunning 
rise in the capabilities of numerical simulation. These 
have shown no indication of local nonlinear rotational 
instabilities (Hawley, Balbus, & Winters 1999), in 
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Keplerian disks. They do, however, reveal nonlinear 
shear instabilities when Corilois forces are absent, or 
when the disk is marginally stable (constant specific 
angular momentum). Indeed, even linear instability is 
possible in some Rayleigh-stable disks, provided that 
global physics is introduced (Goldreich, Goodman, 
& Narayan 1986; Blaes 1987), a result that has been 
numerically confirmed (Hawley 1991). 

The numerical stability findings have been crit- 
icized on the grounds that the effective Reynolds 
number of the codes is too low, and that this damps 
the nonlinear instabilities: the latter require yet-to- 
be resolved spatial scales in order to reveal them- 
selves (Richard & Zahn 1999). In this paper, we show 
that the local disk equations possess a scale invari- 
ance that implies any solution to the governing equa- 
tions must be present on all scales. In other words, 
for every small scale velocity flow, there is an exact 
large scale counterpart with the same long term sta- 
bility behavior. The absence of instability at large 
scale therefore implies the absence at small scales 
as well. Conversely, any true small scale instabili- 
ties (those present in a shear layer, for example), 
must also have large scale counterparts, and there- 
fore instability should be found even at crude numer- 
ical resolutions. This is indeed the case. Our findings 
suggest that if nonlinear hydrodynamical instabilities 
were present in Keplerian disks, such unstable dis- 
turbances must involve dynamics beyond the local 
approximation, and are not an inevitable nonlinear 
outcome of differential rotation. 



2. The Local Approximation 

In cylindrical coordinates {R,(f>,z), the fundamental 
equations of motion for a flow in which viscous effects 
are negligible are mass conservation 



| + v.(H = o, 

and the dynamical equations. 
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Our notation is standard: v is the velocity field, p 
the mass density, P the gas pressure, and $ is the 
Newtonian point mass potential for central mass M: 
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G is the gravitational constant. 

The local limit consists of the following series of 
approximations. First, we assume that R is large and 
z ^ R, so that 
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Choose a fiducial value of R, say Rq. Denote the an- 
gular velocity as fl{R) (we assume a dependence only 
upon i?), and let fio = ^{Ro)- We next erect local 
Cartesian coordinates 

x = R-Ro, y = Ro{^-not) (8) 

which corotate with the disk a.t R = Rq. Let 

w = V — Rilotefj, (9) 

be the velocity relative to uniform rotation ai — 
VIq. In the local approximation, the magnitude of w 
is assumed to be small compared with RQq- 

The undisturbed angular velocity is Keplerian, 



GM 



(10) 



Substitution of equations (|6I10|I into equations (j2HSl 
and retaining leading order, yields the so-called local 
or iJzZZ equations (e.g., Balbus & Hawley 1998): 



|^ + V.(p«;) = 0, 
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The "0" subscript has been dropped in the 2fi terms 

(4) 

in equations 112() and p3|l , and in the derivative term 
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on the right of equation (|12|) . The time derivative is 
taken in the corotating frame, viz.: 

d d d 

— (inertial) = — (corotating) + ^^o^ (15) 

Equations (|11H14|I are well known, and have been 
used extensively in both numerical and analytical 
studies. The fundamental approach dates from nine- 
teenth century treatments of the Earth-moon system 
(Hill 1878). 

3. Scale symmetry in the Hill Equations 

The local equations of motion incorporate an impor- 
tant symmetry in their structure. Let 

wir,t), p{r,t), Pir,t), (16) 

where r = {x,y, z), be an exact solution to the Hill 
equations piH14(l . Then, if a is an arbitrary constant, 

{l/a)w{ar,t), p{ar,t), {l/a^)P{ar,t) (17) 

is also an exact solution to the same equations. The 
proof is a simple matter of direct substitution. 

An equivalent formulation of the scaling symme- 



try is 

w{r/e, t) ^ {l/e)w{r, t) (18) 

p(r/e,t)^p(r,i) (19) 

P{r/e,t) ^ {l/t^)P{r,t). (20) 



In this form, with e <C 1, we see that any solution 
of the Hill equations that involves very small length 
scales has a rescaled counterpart solution with ex- 
actly the same time dependence. In particular, any 
solution corresponding to a breakdown into turbu- 
lence must be present on both large and small scales. 

The implications of this scaling symmetry are of 
particular importance for understanding and test- 
ing the possible existence of local nonlinear instabil- 
ities in Keplerian disks. The key point is that any 
such instability would have to exist not just at small 
scales, but at all scales. Finite difference numerical 
codes would find such instabilities, if they existed. 
Indeed, a constant specific angular momentum pro- 
file is nonlinearly unstable, and is found to be so 
even at resolutions as crude as 32'^. By way of con- 
trast, local Keplerian profiles show no evidence of 
nonlinear instability at resolutions up to 256'^, in- 
stead converge to the same stable solution in codes 



with completely different numerical diffusion proper- 
ties (Hawley, Balbus, & Winters 1999). The argument 
that small scale flow structure is somehow being re- 
pressed is simply untenable. 

To see how the Reynolds number changes with 
scale, assume that a flow is characterized by an ef- 
fective kinematic viscosity v. The scaling argument 
we have just given applies to inviscid equations, so 
we should not expect it to hold in the presence of 
viscosity. The Reynolds number associated with the 
small scale solution is 

w X el , , 

Re, = (21) 

V 

The Reynolds number associated with the large scale 
solution is 



where w here means w{l,t), the value of the veloc- 
ity function evaluated at a fiducial value length / and 
time t. Rei — Rcs/e^ ^ Re, because at larger scales 
both the velocity and the length scales increase by a 
factor of 1/e. In a numerical simulation, strict scal- 
ing invariance is not obeyed. Instead, the large scale 
solutions approach the inviscid limit, while their suffi- 
ciently small scale counterparts are damped. But by 
behaving nearly inviscidly, the large scale solutions 
capture the behavior of the Hill system at all scales. 

4. What this Result Does Not Show 

Obviously, scale invariance does not constitute a 
proof of nonlinear stability in any Keplerian flow. 
There are several points we have not covered. 

First, the local approximation ignores boundary 
conditions. In laboratory flows, the fluid is always 
bounded by hard walls, and boundary layers form. A 
recent laboratory conflrmation of the MRI also finds 
finite amplitude velocity fluctuations in a magneti- 
cally stable flow, for example. But the source of such 
disturbances are boundary layers (Sisan et al. 2004). 

The Hill equations emerge in the limit i? — > oo, 
and therefore curvature terms drop out of the anal- 
ysis. Instabilities that depend, for example, upon 
inflection points or vorticity maxima in the back- 
ground rotation profile would not appear in this limit. 
Nothing precludes them from forming in the w ve- 
locity profile, however, and if such instabilities were 
present they should manifest on large scales as well 
as small. In any case, the criticism of the numerical 
simulations is that extremely small structure is being 



4 



Balbus: Nonlinear Scale Invariance in Local Disk Flows 



lost, and that high Reynolds number differential ro- 
tation is supposedly intrinsically unstable. It is very 
difficult to see how large scale curvature could play an 
essential destabilizing role here. In these equations, 
the curvature terms are nonsingular perturbations. 
Planar Couette and Poiseuille flows break down into 
turbulence without assistance from geometrical cur- 
vature. 

Our Hill analysis together with numerical simu- 
lations would also suggest that a non-Keplerian disk 
with, say, (x R'^-^ is nonlinearly stable. But an 
annulus supporting such a profile is in fact linearly 
unstable (Goldreich, Goodman, & Narayan 1986), 
transporting angular momentum outward even in its 
linear phase. The point is that the annulus supports 
edge modes that become unstable, and these global 
modes do not exist in the local approximation. The 
existence of a similar instabilities in disks found in 
nature cannot be ruled out, though to date none af- 
flicting Keplerian disks have been found. 

The disk thermal structure could also be unstable, 
at least in principle. Nothing presented in this work 
bears on these types of instabilities. 

Finally, there are technical loopholes to the argu- 
ment presented in this paper. What if the unstable 
solution required not just some small scales to be re- 
solved, but very disparate scales? Why this should 
be so is far from clear, but this possibility cannot be 
ruled out a priori. Indeed, one could imagine that 
a fractal structure is required down to inflnitesimal 
scales. Rescaling would not bring such a solution to 
larger characteristic length scales, by definition. This 
solution is obviously not characterized by a critical 
Reynolds number, above which it is necessary to be 
seen. The critical Reynolds number would be infin- 
ity! This is not the argument made by proponents of 
nonlinear high Reynolds number instability. Such a 
solution may remain a mathematical possibility, but 
not one that can be realized in nature. 

5. Conclusion 

The local dynamics of Keplerian or other astrophys- 
ical disk profiles can be captured by a an established 
formalism known as the local, or Hill, approxima- 
tion. The resulting system of equations has an exact 
scale invariance, so that any flow characterized by 
very small scales has an exact large scale counter- 
part with same stability properties. This feature of 
the Hill equations implies that finite difference codes 
at available resolutions are sufficient to explore the 



possibility of local nonlinear shear instabilities in as- 
trophysical disks. If simulations accurately describe 
the large scale behavior of the Hill system, there is 
nothing more to uncover at small scales; it is simply 
renormalized large scale behavior. The absence of any 
observed instabilities in Keplerian numerical studies, 
coupled with the ready manifestation of such insta- 
bilities in local shear layers and constant specific an- 
gular momentum systems, suggests that any putative 
nonmagnetic disk instability would have to incorpo- 
rate physics beyond simple differential rotation. 
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